Loading the data

The file Market.rda contains the list Market.

The first element of the list Market is an xts object corresponding to the value of the SP500 index from “2000-01-03” until “2013-09-10”.

The second element vix corresponds to the value of the VIX index at the same dates.

The third element rf corresponds to the term structure on “2013-09-10” for 14 maturities varying from 1 day to 30 years.

The fourth (calls) and the fifth (puts) contains the strikes (K), maturities (tau) and implied volatilities (IV) for calls and puts options respectively.

We plot the volatility surface in Figure .

Volatility surface for the put options in the list Market.

Pricing of a portfolio of options

The prices of the call options using the Black-Scholes formula are given in table \ref{tab: calls_prices)

Call options prices
S0 K tau sigma r C
1683.99 1600 0.08 0.1453 0.00115 87.569
1683.99 1650 0.08 0.1453 0.00115 47.724
1683.99 1750 0.16 0.1453 0.00144 15.302
1683.99 1800 0.16 0.1453 0.00144 6.380

The price of this book (the sum of prices of the calls that constitutes the book) is \(V_0\) = 156.98.

One risk driver and Gaussian model

In this section, we assume that the risk of our book of options is driven by a single factor namely the change in price of the SP500 index. The logreturns of the index are assumed to follow a normal distribution with parameters given by the mean and standard deviation from the historical values.

\label{fig:PLdist_gauss1D} P&L distribution for normal log returns of the SP500. The red line represents the value at risk at 95% probability.

P&L distribution for normal log returns of the SP500. The red line represents the value at risk at 95% probability.

We model 10000 possible values for the stocks in 5 days and use them to estimate the distribution for the P&L in a week (5 days). The distribution of the P&L is given in Figure .

The VaR at risk at the 95% confidence level is VaR95 = 122.8 and the expected shortfall is ES95 = 134.97. So we expect the book of options to lose at least 122.8 from its initial value 5% of the time.

Two risk drivers and Gaussian model

In this section, we assume that the risk of our book of options is driven by two factors namely the change in price of the SP500 index and the volatility of the market (given by the vix). The logreturns of the SP500 and the vix are assumed to follow a normal distribution with parameters given by the mean and standard deviation from the historical values.

\label{fig:PLdist_gauss2D} P&L distribution for normal log returns of the SP500 and the vix. The red line represents the value at risk at 95% probability.

P&L distribution for normal log returns of the SP500 and the vix. The red line represents the value at risk at 95% probability.

We model 10000 possible values for the SP500 and the vix in 5 days and use them to estimate the distribution for the P&L in a week (5 days). The distribution of the P&L is given in Figure .

The VaR at risk at the 95% confidence level is VaR95 = 111.63 and the expected shortfall is ES95 = 123.73. So we expect the book of options to lose at least 111.63 from its initial value 5% of the time.

The results don’t differ much from the case where only the stock was considered to drive the risk of the book.

Two risk drivers and copula-marginal model (Student-t and Gaussian copula)

In this section, we assume that the risk of our book of options is driven by two factors namely the change in price of the SP500 index and the volatility of the market (given by the vix). The logreturns of the SP500 are assumed to follow a Student-t distribution with \(\nu = 10\) degrees of freedom and the logreturns of the vix are assumed to follow a Student-t distribution with \(\nu = 5\) degrees of freedom. A normal Copula is assumed to merge the marginals.

\label{fig:PLdist_ststgauss2D} P&L distribution for normal log returns of the SP500 and the vix. The red line represents the value at risk at 95% probability.

P&L distribution for normal log returns of the SP500 and the vix. The red line represents the value at risk at 95% probability.

We model 10000 possible values for the SP500 and the vix in 5 days and use them to estimate the distribution for the P&L in a week (5 days). The distribution of the P&L is given in Figure .

The VaR at risk at the 95% confidence level is VaR95 = 99.07 and the expected shortfall is ES95 = 113.1. So we expect the book of options to lose at least 99.07 from its initial value 5% of the time.

Volatility surface

We model 10000 possible values for the SP500 and the vix in 5 days and use them to estimate the distribution for the P&L in a week (5 days). The simulated values of the vix are then shifted by the difference between the one year ATM value of the implied volatility and today value of the vix. The distribution of the P&L is given in Figure .

The VaR at risk at the 95% confidence level is VaR95 = 108.64 and the expected shortfall is ES95 = 121.19. So we expect the book of options to lose at least 108.64 from its initial value 5% of the time.

Full approach

Here, we fit a GARCH(1,1) to the log returns of the SP500 and a AR(1) process to the log returns of the VIX. The residuals of the two processes are assumed to be invariants and to be linked by a Gaussian copula. We generate 10000 possible values for the SP500 and the vix in 5 days and use them to estimate the distribution for the P&L in a week (5 days). The distribution of the P&L is given in Figure .

The VaR at risk at the 95% confidence level is VaR95 = 14.98 and the expected shortfall is ES95 = 16.32. So we expect the book of options to lose at least 14.98 from its initial value 5% of the time.